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Abstract. We present the first spatially resolved spectroscopic observations of the recently discovered quadruple QSO and 
gravitational lens HE 0435 — 1223. Using the Potsdam Multi-Aperture Spectrophotometer (PMAS), we show that all four QSO 
components have very similar but not identical spectra. In particular, the spectral slopes of components A, B, and D are in- 
distinguishable, implying that extinction due to dust plays no major role in the tensing galaxy. While also the emission line 
profiles are identical within the error bars, as expected from lensing, the equivalent widths show significant differences between 
components. Most likely, microlensing is responsible for this phenomenon. This is also consistent with the fact that component 
D, which shows the highest relative continuum level, has brightened by 0.07 mag since Dec 2001. We find that the emission 
line flux ratios between the components are in better agreement with simple lens models than broad band or continuum mea- 
surements, but that the discrepancies between model and data are still unacceptably large. Finally, we present a detection of the 
lensing galaxy, although this is close to the limits of the data. Comparing with a model galaxy spectrum, we obtain a redshift 
estimate of zjcns = 0.44 ± 0.02. 
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1. Introduction 

HE 0435—1223 was first discovered in the course of the 
Hamburg/ESO survey for bright QSOs dWisotzki et alj EoOO). 
and recently found to be a rather spec tacular example of 
a quadruply imaged QSO dWisotzki et aL .2002; hereafter 
Paper I). The quasar is at a redshift of z = 1.689 and had 
a total magnitude g — 17.8 at the epoch of discovery, with 
evidence for significant variability. The conspicuous geometry 
of four blue point sources of similar brightness arranged with 
nearly perfect symmetry around a red elliptical galaxy left no 
doubt about the nature of this object as being gravitationally 
lensed. However, Paper I featured only a blended spectrum of 
all components, and a formal verification based on individual 
spectroscopy of each component is still pending. In this pa- 
per we present the first set of spatially resolved spectra of this 
quadruple QSO, obtained by integral field spectrophotometry. 

2. Observations and data reduction 

We targeted HE 0435—1223 during the first regular (i.e. non- 
comissioning) observing run of the Potsdam Multi-Aperture 
Spectrophotometer PMAS, mounted at the Calar-Alto 3.5 m 
telescope for several nights between 2-7 September, 2002. 
Detai ls of this new instrument are given elsewhere ( Rot h et alJ 
I2OO0I) : here we only summarise its capabilities. PMAS cur- 
rently consists of a 16 X 16 elements microlens array (an up- 



grade to 32 X 32 elements is planned) coupled by optical fi- 
bres to a purpose-designed spectrograph. The image scale is 
0'.'5 per spatial pixel {'spaxel'), thus the total field of view is 
8" X 8". We used the V300 grating, giving a spectral reso- 
lution of 6A FWHM and a spectral range of 3300 A, 
from 3950 to 7250 A in our preferred setting. The data were 
recorded by a SITe 2048 x 4096 pixels CCD detector, read 
out in 2 X 2 binned mode. PMAS also features an independent 
on-axis imaging camera with a 3.'4 x 3.'4 field of view and a 
1024 X 1024 pixels detector, mainly used for acquisition and 
guiding (termed the AG camera henceforth). 

At the observing date in early September, HE 0435 — 1223 
became visible only very shortly before morning twilight, and 
only at airmasses of around 2. We revisited the object repeat- 
edly and could record useful data at the end of three differ- 
ent nights, with exposure times of altogether 3 x 1800 s and 
1 X 1200 s. The external (DIMM) seeing was usually around or 
slightly below 1", which translated into an effective seeing in 
our PMAS data of between I'.'l and 1'.'3. Wavelength calibra- 
tion and tracing of the 256 fibre spectra were facilitated through 
exposures using an internal calibration unit before or after each 
science spectrum. 

Observations of the spectrophotometric standard star Hz 4 
in the morning twilight provided the data needed for flux cal- 
ibration. Since the AG camera was equipped with a Johnson 
V filter, we could also exploit the acquisition exposures to 
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Fig. 1. Top row: Example PMAS data 
of HE 0435-1223, for different wave- 
bands. The first three images are quasi- 
monochromatic (3.3 A bandwidth), the 
right-hand image is averaged over a 
wide passband. Each image measures 
8" X 8"m and each spatial pixel is 
O'.'S X 0'.'5. The orientation is stan- 
dard (North is to the top. East is 
to the left). The four QSO compo- 
nents are clearly resolved. Notice the 
strong shift of image centroids as a 
function of wavelength (cf. Fig. |2ll, 
caused by differential atmospheric re- 
fraction. Middle row: CoiTesponding 
model images consisting of four super- 
posed Gaussians (for details see text). 
Bottom row: Residuals after subtract- 
ing the four point sources. 



check and improve the spectrophotometric calibration accu- 
racy. Despite the high airmass of HE 0435 — 1223 during the 
observations, the flux level measurements of corresponding 
quasar spectra taken in different nights were remarkably con- 
sistent, with deviations of the order of 1 %. This indicates both 
photometric observing conditions and a stable system. Notice 
that since the PMAS microlens array reimages the exit pupil of 
the telescope, there are no geometrical losses due to incomplete 
filling of the focal plane (as with some other integral-field in- 
struments). PMAS is thus certainly suitable for accurate spec- 
trophotometry, and we shall use this capability below. 

The data were reduced u sing our own IDL-based software 
package P3d ('Becker''2002). The reduction consists of stan- 
dard steps such as debiasing and flatfielding using twilight ex- 
posures, and dedicated routines such as tracing and extracting 
the spectra of individual fibres and reassembling the data in 
form of a three-dimensional data cube. Wavelength calibration 
and rebinning to a constant spectral increment is also part of the 
reduction procedure. The top row of Fig.^shows as examples 
some quasi-monochromatic and synthetic broad-band images, 
giving a good representation of the image quality, in particular 
S/N level and angular resolution. In fact, the S/N measured in 
our PMAS data is fully consistent with expectations from pure 
shot noise. 



3. Deblending of the QSO components 

3. 1. Procedure 

Starting with the data cubes produced on output by the P3d 
package, we developed a simple but powerful scheme to extract 
all four QSO spectra simultaneously. The following considera- 
tions guided our sti-ategy: 



- Because of the small field of view, the point-spread function 
was not known a-priori, but had to be determined within 
the QSO data.' Non-simultaneous observations of PSF cal- 
ibrator stars were not considered an option, because of the 
significant temporal PSF variability at the site. 

- Due to differential atmospheric refraction, the source cen- 
troids vary with wavelength relative to the grid of 16 x 16 
spatial pixels. This effect is particularly grave for our high- 
airmass data, leading to a shift of more than four pixels (i.e. 
> 2") over the spectral range covered (see Figs.^andlJ}- 

- Over the substantial spectral range covered, the width of 
the PSF can not assumed to be constant. This effect turned 
out to be less relevant than expected in our data, cf. Fig. 
|3l but significantly affected other data taken in the same 
campaign. 

Our adopted algorithm is conceptually similar, albeit somewhat 
more complex, to procedures already used by us and others for 
the simult aneous extraction o f double QSO spectra in long-slit 
data (e.g..lSmetteetal.ll995h. 

In essence, the routine performs a multi-component fit of 
n two-dimensional Gaussians plus a spatially constant back- 
ground value to the data in each quasi-monochromatic image, 
using a modified Levenberg-Marquard minimisation scheme 
JPress et al. 1992). In the case of HE 0435-1223, we adopted 
n = 4 for the four QSO components. At the chosen spectral 
range, it was not strictly necessary to include also the lensing 
galaxy; as discussed in Sect. Qbelow, we have marginally de- 
tected the galaxy, although on a level not much above the PSF 
mismatch residuals. 

' In the future it will be possible to use the guide star data recorded 
by the AG camera during a PMAS exposure to establish a simultane- 
ous PSF reference, but this mode was not yet operational at the time 
of observing. 
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Fig. 2. Centroid displacement as a function of wavelength, 
mainly due to differential atmospheric refraction. Each point 
corresponds to one quasi-monochromatic spectral layer. In the 
later fitting passes, the measured centroids are replaced by the 
polynomial approximation plotted by the solid lines. 



On output, the procedure provides the fitted parameters and 
their standard errors in a table. Each Gaussian is basically char- 
acterised by six parameters per spectral element; centroids in x 
and y, FWHM along major and minor axis, position angle, and 
amplitude. By including some prior knowledge we can reduce 
the number of free parameters: First, the values of FWHM are 
assumed to be the same for all four point sources. The pro- 
cess is then iterated: In an initial pass, all parameters are fitted 
freely. In subsequent passes, the centroids, FWHM values and 
position angles are replaced by smoothly varying polynomial 
functions of A, as exemplified by Figs. Inland [S] Once fitted, 
these parameters are held fixed and do not participate in the fit- 
ting process, leaving only the 4 amplitudes to be fitted in the 
final pass. Model and residual data cubes are constructed along 
for inspection (see Fig.[3. The procedure turned out to be ro- 
bust and efficient, simultaneously providing deblended spectra 
and a reasonable PSF, as demonstrated by the small residuals 
in the bottom row of Fig.^ 

After the final pass, the product of fitted amplitudes and 
(FWHM)^ gave the extracted number of counts per spectral 
element for each QSO component. A similar extraction per- 
formed for the spectrophotometric standard star Hz 4 provided 
a formal flux calibration. We note that the efficiency of PMAS 
at 4000 A is stin 30 % of its peak value at 6000-7000 A, which 
is remarkably high for an integral field spectrograph. 

3.2. Results 

Visual inspection of the images confirms that the four com- 
ponents are well separated, with little or no significant mu- 
tual contamination (see Fig.[l). The extracted spectra from the 
four exposures of HE 0435—1223, obtained during three differ- 
ent nights, were combined using an inverse variance weighting 
scheme, which was also used to eliminate bad data values con- 
taminated by cosmics or other spurious signals. Figurel^shows 
the four final spectra together with the associated error vectors. 
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Fig. 3. Same as Fig.|2l but for the full width at half maximum 
(average of major and minor axis) of the Gaussian PSF. 

The S/N in the continuum around A = 5000 A is about 60 for 
component A and ^ 40 for the three fainter components B, C, 
andD. 

Apparently, all four components show very similar and typ- 
ical quasar spectra. While this is not really a surprise given 
the configuration of the system and the colours measured in 
Paper 1, it is nevertheless the first time that individual spec- 
tra of the components of HE 0435 — 1223 are published. The 
redshifts derived from the C IV and C III] lines are identical to 
within the measurement accuracy: z — 1.6895 ± 0.0005; the 
redshift scatter between the components is Az ~ 0.0015 (rms), 
corresponding to less than one spectral pixel. 

4. Spectrophotometric Analysis 

4. 1. Broad-band fluxes 

Our broad-band data from Paper 1 were taken with the 
Magellan telescope in SDSS filters g, r, and i at two epochs, 
14 Dec 2001 and 12 Feb 2002. In order to facilitate a sti-aight- 
forward comparison between these and our new PMAS data, 
we list synthetic broad-band magnitudes in Table [fl obtained 
by integrating over the g and r filter curves (our spectra do 
not cover the i band). At the high S/N of our data, the for- 
mal errors of the synthesised broad-band values are far below 
0.1 % and presumably outweighed by unknown systematics. 
The overall flux level in the g band increased by 0.5 mag be- 
tween Dec 2001 and Sep 2002. Already in Paper 1 we noted 
that HE 0435 — 1223 was variable, having brightened by almost 
0.2 mag between Dec 2001 and Feb 2002. We did not detect 
variability on timescales of days: The fluxes measured in the 
three consecutive observing nights during this run are constant 
to within ^ 1 %. 

The magnitude differences between the components are 
largely consistent with the values quoted in Paper 1, with one 
exception: Component D has brightened by 0.07 mag, from be- 
ing the faintest component in the Magellan data to now at the 
same level as component B. We argue below that this is most 
likely due to microlensing. 
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Fig. 4. Extracted and coadded spectra of the four QSO components. The thin Unes near zero represent the Icr eiTor arrays. The 
adopted local continuum levels, used to derive emission line properties, are indicated by the dashed line segments. 



We can also provide the first reliable colour measurements. 
The colours listed in Paper 1 were based on a very crude cal- 
ibration scheme, tied to just one star with rather inaccurately 
known intrinsic colours. We find that the colours for all four 
components are basically identical, {g — r)ABCD = 0.35 with 
an rms scatter of 0.02 mag. In Paper 1 we found component C 
to be just slightly redder than the mean by 0.03 mag; although 
this is clearly at the edge of our measurement accuracy, we find 
that {g — r)c = 0.38, and the colour deviation is exactly repro- 
duced in our new spectra. However, these results do not nec- 
essarily imply differences in continuum slope, as broad-band 
colours may well be influenced, at this accuracy level, by the 
different emission line equivalent widths discussed in the next 
section. 

4.2. Emission lines 

An important spectrophotometric diagnostic for gravitationally 
lensed QSOs is the comparison of emission line properties be- 
tween components. While ideally the spectra of lensed quasar 
components should look identical, microlensing and other con- 
taminating effects may lead to differences. In order to derive 
the emission line properties, we had to define the underlying 
continuum. Since we were not interested in a global continuum 
level, we defined two ad hoc continuum windows left and right 



of each line and connected these by straight lines (see Fig.|4}. 
The windows were located at [3980 A, 4020 A] and [4530 A, 
4630 A] for the C IV A 1550 line, and at [4890 A, 4950 A] and 
[5380 A, 5440 A] for the Allll/Clll] A1909 blend. The local 
continua thus derived are indicated by the dashed line segments 
in Fig. 0] This simple procedure worked quite well, although 
the continuum level at the blue side of the C IV line is poorly 
constrained due to the short wavelength cutoff of the spectra. 

We then integrated the continuum-subtracted spectra within 
narrow windows centred on the emission line cores. The selec- 
tion of these windows was to some extent arbitrary and mainly 
guided by the wish to avoid the low S/N wings which are most 
strongly affected by errors in the continuum estimate. The fi- 
nally adopted integration windows were [4020 A, 4300 A] for 
Civ and [5075 A, 5 180 A] for Cm]. Table 13 lists the results 
for all four components. We do not quote formal statistical er- 
rors because, as before, these are very small, typically ~ 1 % 
for the line fluxes, and almost certainly dominated by unknown 
systematical errors. 

With high significance, the equivalent widths of the lines 
are not identical in the four components. This is akeady ap- 
parent from visual inspection of the spectra in Fig.E] In other 
words: The flux ratios between the QSO components are dif- 
ferent in the continuum and in the broad emission lines. We 
discuss this phenomenon in the next sections. 
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Table 1. Synthetic broad band magnitudes for HE 0435-1223, 
in the Johnson V and in the SDSS g and r bands (Vega magni- 
tudes). Labelling of components as in Paper 1, i.e. clockwise. 



Component 


V 


9 


r 


g~r 


A 


18.34 


18.50 


18.16 


0.34 


B 


18.88 


19.03 


18.70 


0.33 


C 


18.96 


19.15 


18.77 


0.38 


D 


18.90 


19.07 


18.72 


0.35 


A+B+C+D 


17.23 


17.40 


17.05 


0.35 



Table 2. Emission line fluxes and rest-frame equivalent 
widths for the four components. Fluxes are expressed in 
equivalent widths in A. 



10 ^® erg s ^ cm ^ 





Civ A1550 


Cm] 


A 1909 


Component 


/ 


Wo 


/ 


Wo 


A 


120.1 


19.2 


60.3 


13.5 


B 


90.4 


24.2 


46.3 


17.2 


C 


86.0 


26.8 


42.3 


16.9 


D 


66.2 


18.3 


31.7 


11.8 



Finally, we also compared the emission line profiles. To this 
purpose we rescaled all continuum-subtracted lines to a com- 
mon line integral and superimposed them pairwise. This is doc- 
umented in Fig.|5] where the normalised C IV and C III] lines 
of each of the components B, C, and D are plotted over A. As 
expected for a gravitationally lensed QSO, the individual line 
profiles are very similar and consistent with the assumption that 
they originate in the same object. 

5. Spectral differences: Evidence for 
Microiensing 

In the previous section we documented that despite the certain 
nature of HE 0435—1223 as a lensed QSO, the spectra of indi- 
vidual components are not strictly identical. Several effects can 
be invoked to explain spectral differences between lensed im- 
ages: Foreground extinction, intrinsic variability, or microiens- 
ing. We discuss each of these in turn. 

Extinction is probably not important in this object. The 
nearly identical (and typically quasar-like) colours of the four 
components imply that there is no significant differential ex- 
tinction between these lines of sight. Although it is conceivable 
that all lines of sight suffer from the same amount of extinc- 
tion, this would be an extreme coincidence if the dust column 
density were high. Furthermore, the lensing galaxy is a large 
elliptical (cf. Paper 1), not expected to contain large quantities 
of smoothly distributed dust in its outskirts. And of course the 
different equivalent widths of the emission lines cannot be ex- 
plained by dust extinction at all, so there has to be at least one 
additional effect. 

Intrinsic flux variability is certainly an issue in 
HE 0435 — 1223, as already discussed in Sect. 14. ll However, 
it is unlikely that variability is the explanation for the dif- 
ferences. From the lensing geometry, the light travel time 



delay between the images is predicted to be of the order of 
a few days (Paper 1). This means that very rapid and drastic 
intrinsic spectral changes would be needed in order to create 
differences between the images; we see no evidence at all 
for such violent variability. Even if there were such changes, 
escaping our attention, these would have to appear in each of 
the components consecutively within days. The fact that the 
spectra taken at three different nights look identical within the 
error bars (i.e. with an uncertainty of no more than a few per 
cent) speaks against significant short-time spectral variability. 

Gravitational microiensing remains as the most straightfor- 
ward explanation for the observed spectral differences between 
the four images. Within that scenario, microiensing affects the 
continuum and not the emission lines because of the sizes in- 
volved: The broad-line region in quasars is believed to be much 
larger than the typical cross-section of stellar mass lenses, and 
any magnification patterns are thus completely smeared out. 
On the other hand, the continuum-emitting region is orders of 
magnitude smaller than the BLR and certainly prone to mi- 
croiensing (de)magnifi cation, as most drama ticallv seen in the 
'Einstein C ross' (e.g., IW ozniak et al. 2000). Spectroscopy of 
that quasar jLewis et all 1998.) showed also differences in the 
emission line equivalent widths very similar to the case dis- 
cussed here. The most extreme known case of selective contin- 
uum amplificatio n bv microiensing is t he 'Double Hamburger' 
HE 1104-1805 dWisotzki et alJll993h . where there is even a 
strong chromatic effect, i.e. a wavelength-dependent amplifica- 
tion of the continuum (the spectrum gets bluer because the con- 
tinuum region is smaller at shorter than at longer wavelengths 
and hence is more magnified; cf. IWambsganss & Paczvnskil 

dm). 

For all these objects as well as for HE 0435 — 1223, the 
emission line profiles are identical in all components. The spec- 
tral differences can then be interpreted as being entirely due 
to additional continuum contributions on top of the intrinsic 
'macrolensed' spectra. Since microiensing is time-variable as a 
result of non-negligible effective transversal motions, this addi- 
tional contribution is expected to vary. Interestingly, it is com- 
ponent D for which we find evidence for a relative brightening 
between Dec 2001 and Sep 2002; that component has the low- 
est equivalent widths and thus, in the microiensing picture, the 
highest continuum amplification. We conclude that microiens- 
ing is the only plausible explanation for the observed differ- 
ences in line-to-continuum ratios. On the other hand, there 
is currently no indication of significant differential chromatic 
continuum amplification, since all components have so similar 
colours. 

At this point it is useful to reconsider briefly the question 
whether the emission lines might also be affected by microiens- 
ing. The broad-line region in quasars (BLR) is generaUy held 
to be much bigger than the continuum region; for a quasar of 
the intrinsic (de-lensed) luminosity of HE 0435—1223, a BLR 
radiu s of the order of l OOlight-davs can be expected (Kaspi 
et al. teaspi et al."2000')), two orders of magnitude larger than 
the Einstein radius of a solar-mass star at z ~ 0.4. Any mi- 
croiensing would therefore act on a small region of the BLR 
only, possibly causing some variations of the line profile as an 
observable signature, but certainly no significant overall am- 
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Fig. 5. Pairwise comparison of nor- 
malised emission line profiles (upper 
row: C IV; lower row: C III]). The light 
thick lines represent component A, the 
dark thin lines denote components B, 
C, and D, respectively. Below each 
panel we show also the difference be- 
tween the profiles. 



Table 3. Observed and predicted flux ratios between the four 
components: (1) CiV and Cm] emission lines; (2) contin- 
uum ratios at A = 4580 A; (3) broad-band measurement from 
Paper 1 ; (4) flux ratios predicted from the model in Paper 1 . 



Ratio 


lines 


continuum 


g band 


Model 


B/A 


0.77 


0.62 


0.60 


1.11 


C/A 


0.71 


0.54 


0.57 


1.04 


D/A 


0.56 


0.61 


0.55 


0.68 


C/B 


0.93 


0.87 


0.95 


0.94 


D/B 


0.73 


0.98 


0.93 


0.61 


D/C 


0.79 


1.12 


0.97 


0.65 



plification (cf. also lSchneider & Wambsganssll990l) . Even that 
effect cannot be strong in HE 0435 — 1223, since the line pro- 
files are so similar (cf. Fig.|5}. We therefore assume in the fol- 
lowing that a clean distinction between the two domains exists: 
While the continuum source is very likely significantly ampli- 
fied by microlensing, the emission lines remain more or less 
unaffected. 

6. Flux ratios 

One interesting consequence of the microlensing interpretation 
is the expectation that the emission line fluxes should mirror 
the 'true' flux ratios of the four components, i.e. those given by 
macrolensing of the smooth galaxy potential alone. This is im- 
portant for the construction of mass models for the lens: While 



flux ratios from broad-band photometry are notoriously unreli- 
able because of the ever present possibility of microlensing, we 
suggest that emission line flux ratios can be directly used as ad- 
ditional model constraints. In Table|3]we provide the pairwise 
emission line flux ratios between the components. As a kind of 
'sanity check', we also determined the corresponding contin- 
uum flux ratios, listed in the second column of Table |3] Since 
the spectral slopes are so similar, these values do not change 
significantly over the spectral range covered. 

In Paper 1 it was found that a simple model involving a 
singular isothermal sphere plus external shear could reproduce 
the positions of the QSO components and of the lensing galaxy 
well within the error bars. However, the model (which was 
based on astrometric constraints only) predicted components 
B and C to be substantially brighter, by more than 0.5 mag 
relative to A and D, than observed in broad-band colours (cf. 
the last two columns in Table |3}. Such disagreement between 
models and data has been seen in several other lensed QSOs 
as well, and two hypotheses have recently been proposed to 
account for these anomalies. On the one hand, they may indi- 
cate the presence of su bstructure in the lensing potential (e.g., 
iMetcalf &Zhaoll2002h : or they could be due to microlensing 
at large optical depths, especially by suppressing the flux of 
saddle point images (iSchechter & Wamb s_ g_anss 2002). These 
two explanations differ significantly in their prediction of the 
flux ratios of extended source components. While the flux ra- 
tios should not depend strongly on source size in the case of 
a lens with substructure, a necessary consequence of the mi- 
crolensing hypothesis is that the flux ratios should approach 
the 'macrolens' values as soon as the source is too extended 
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to be affected by the microlens caustic network. This impHes 
that if the size of the QSO broad-line region is already in this 
domain, emission-line flux ratios should be much closer to the 
model values if the model is reaUstic. There is at least one 
known quadruple QSO where this se ems to be the case, namely 
B 1422+231 (cf. llmDev eta im^). 

This test can be performed with our spectrophotometric 
measurements of HE 0435 — 1 223 . Table|3]shows that the agree- 
ment between measured flux ratios and the model predictions 
is indeed much better if one uses the emission line values than 
with either broad band or continuum data. Components B and 
C are brighter in the emission lines than in the continuum, rel- 
ative to A and D, while the ratios A/D and B/C do not change 
significantly. These changes are precisely in the direction ap- 
proaching the model, just as predicted in the microlensing sce- 
nario by Schechter & Wambsganss. 

However, the agreement is anything but good, and espe- 
cially component A is still much brighter than is expected. 
As this could be due to shortcomings of the simple model 
which did not even take flux ratio constraints into account, 
we attempted to compute a new model. We used C. Keeton's 
software tool lensmodel lKeeton.2001 ) to explore the con- 
sequences of adding measured flux ratios as additional con- 
straints. For consistency, we stayed within the singular isother- 
mal sphere plus external shear model (although we also tried 
some elliptical models). To our surprise, it was extremely dif- 
ficult to make the predicted flux ratios approach the measured 
ones; the predictions tended to stay close to the Paper 1 values, 
which were based on the astrometry of the four QSO compo- 
nents alone (the position of the lensing galaxy was not con- 
sidered to be sufficiently constrained, and therefore it was al- 
ways left to float, as in Paper 1). This would only change when 
the flux ratio constraints were taken to be much tighter than 
the astrometric constraints. For example, we could roughly fit 
the flux ratios when allowing for several tenths of arc sec- 
onds of positional shifts, corresponding to a formal deviation 
of the order of more than ~ SOcr, which is totally unacceptable 
given the excellent astrometric quality of the Magellan data. 
Thus, even using the emission line measurements we failed to 
build a consistent model for the lensing potential that correctly 
reproduces both positions and flux ratios. But we recognise 
that a more sophisticated model of the lens mass distribution 
might well prove this statement premature. We have deliber- 
ately limited ourselves to explore only a narrow range of plau- 
sible configurations and now invite the experts in the field to 
follow up on us. Until then, we conclude that our measure- 
ments give some support to the microlensing hypothesis by 
ISchechter & W ambsganss ('2002^, but that no wholly satisfy- 
ing explanation for the flux ratios in HE 0435 — 1223 can yet be 
given. 

7. The lensing galaxy 

The Magellan i band image presented in Paper I shows the 
lensing galaxy prominently, close to the geometric centre be- 
tween the four QSO images. A similarly clean detection in 
our new PMAS data was not to be expected, for two reasons: 
(i) The much coarser pixel grid (0'.'5 for PMAS vs. 0'.'07 for 
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Fig. 6. Left panel: 4x4 supersampled image of 
HE 0435-1223, coadded over a simulated 5500-7200 A 
passband (see text for details). Right panel: Residual image 
after PSF subtraction, supersampled and coadded in the same 
way. Despite the significant PSF residuals, the lensing galaxy 
is clearly detected. The positions of the four QSO components 
and the lensing galaxy (from Paper 1) are marked by crosses. 
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Fig. 7. Redshift estimation for the lensing galaxy. Left panel: 
vs. redshift for model fit of 10 Gyr SSP to the binned resid- 
ual spectrum at position of the lensing galaxy. Right panel: 
Residual spectrum with best-fit SSP model spectrum superim- 
posed. For details see text. 



MagellanH-MagIC); (ii) the inferior effective seeing I'.'l 
here vs. 0'.'6 for Paper I). Furthermore, the nonavailability of a 
simultaneously obtained PSF reference made a clean PSF sub- 
traction difficult. We have nevertheless attempted to detect the 
lensing galaxy, both in 'imaging' and in 'spectroscopy' mode. 
In our search we made use of some pieces of information avail- 
able from Paper I. Firstly, we know the location of the galaxy 
relative to the QSO point sources. Secondly, we assume that its 
redshift is located within 0.3 ^ z ^ 0.5, to be in accordance 
with the photometric redshift estimate in Paper I. Notice that 
the narrowing of the range to0.3^z^0.4in Paper I using a 
luminosity-mass argument should be disregarded, as there was 
an error in computing the rest-frame absolute magnitude. 

In order to partly overcome the coarse pixel grid limitation, 
we have exploited the fact that due to atmospheric refraction, 
the QSO centroid gradually shifts as a function of wavelength. 
In a way, this can be seen as almost equivalent to extensive 
dithering with respect to a broad-band image. The image is 
shifted by several pixels over the entire wavelength range, a 
fortuitous byproduct of the high airmass during observation. 
Within a band of just ~ 1000 A, the displacement is already 
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a good fraction of a pixel. Furthermore the number of samples 
taken is huge, roug hly 300 for a 1000 A bandpass. The infor- 
mation volume required for substantial subsampling therefore 
certainly exists. We have written a small application to coadd 
monochromatic data planes into a broad-band image, including 
a 4 X 4 subsampling and applying the polynomial approxima- 
tion to the centroid shifts described in Sect. 13. II 

Guided by the redshift prior on the lensing galaxy, we se- 
lected only the red part of the spectrum, 5500 A < A < 7200 A. 
In addition, pieces of the spectrum heavily contaminated by 
night sky lines were left out. Both the original and the residual 
data cubes were run through this procedure. The result is shown 
in Fig. |6l The left panel shows that the four point sources are 
now much better defined than in the original coarse pixel data. 
In the right panel we find significant PSF residuals, featuring 
indications of positive central elevations and negative anullus- 
like patterns which document the limits of our simple assump- 
tion to model the PSF as a Gaussian. But we also find that the 
strongest residual is always positive and occurs very close to 
the expected position of the lensing galaxy. Relative to the po- 
sition of component A, the difference from the value quoted 
in Paper 1 is only ~ O'.'l in right ascension and 0'.'2 in decli- 
nation. Going through the same exercise with a A < 5500 A 
image does not yield any discernible signal at that location. We 
therefore conclude that the detection is real and significant. 

Thus encouraged, we looked whether spectral information, 
possibly even the lens redshift could be derived. We there- 
fore extracted and coadded the spectra of a 2 x 2 pixel block 
(1" X 1") at the expected location from each residual data cube. 
The result was an extremely noisy spectrum, with a signal that 
was on average positive above 6000 A (S/N of ^ 2, for the orig- 
inal 3.3 A wavelength steps) and zero or slightly negative be- 
low that. No QSO emission lines were present, indicating that 
cross-talk from the QSO point sources was negligible. Since it 
was hopeless to detect stellar absorption lines from the lensing 
galaxy in such low S/N data, we degraded the spectral resolu- 
tion by summing the data into bins of 300 A width, with the in- 
tention to use the overall spectral energy distribution for a red- 
shift estimate. This is documented in Fig.0 We assumed an old 
stellar population to dominate the SED of the lens (which has 
elliptical morphology, cf. Paper I), and cross-correlated the data 
with a solar m etallicity, age of 10 Gyrs model galaxy spectrum 
from jjimenez et al. 2000) . kindly provided by R. Jimenez. 

The left panel of Fig. [T] shows that x^(^) has a strong but 
flat-bottomed minimum between z — 0.43 and z — 0.45, 
and a turnover at z = 0.58 towards another local minimum 
at much higher z. A lens redshift z > 0.5 is not compati- 
ble with our prior assumptions; apart from its incompatibility 
with the colour measurement of Paper I, such a redshift would 
make the lens extremely massive. We therefore settle at a red- 
shift estimate of zicns = 0.44 ± 0.02, for which we show our 
best fit model overplotted over the binned residual spectrum. 
Note that we obtain consistent results (within the quoted error 
range) when analysing the individual exposures taken in differ- 
ent nights. 

The model roughly reproduces the main feature of the data, 
its steep rising towards the red, but otherwise the fit is not good, 
and the error bars are probably underestimated. Selecting a 



younger stellar population would further degrade the fit quality, 
since already now our data are too red compared to the 10 Gyrs 
model (which is already the age of the universe at z = 0.44). 
We conclude that our new lens redshift is a significant improve- 
ment over the rough colour estimate of Paper I, but the value 
is nevertheless tentative. A measurement of zicns involving the 
detection of stellar absorption lines in the lensing galaxy is still 
indispensable. 

8. Conclusions 

We have presented the first spatially resolved spectra of the new 
quadruple QSO HE 0435 — 1223. The choice of an integral field 
spectrograph, and in particular the design of PMAS, enabled us 
to obtain high signal-to-noise data of spectrophotometric qual- 
ity for all four quasar images simultaneously. The increasing 
availability of such instruments at many major telescopes opens 
a new window of opportunity to very efficiently observe grav- 
itationally lensed QSOs - especially quadruple systems which 
are very hard to observe with traditional slit spectroscopy. Yet, 
two limitations need to be overcome to make integral field sys- 
tems fully superior to other spectrographs, as far as gravitation- 
ally lensed quasars are concerned. 

Firstly, resolving close systems requires an accurate knowl- 
edge of the point spread function. This is easily supplied by 
modern imaging cameras and multi-slit spectrographs, but the 
field of view of existing integral field units is generally much 
too small to simultaneously capture an isolated PSF reference 
star. One possible way out would be to use a second imag- 
ing device in parallel, such as is planned with the AG camera 
of PMAS. Another option might be to implement clever 'PSF 
self-calibration' algorithms, alleviating the need for an external 
PSF reference. Our simple approach involving just an analytic 
expression could be no more than a preparatory study to such 
techniques. 

Secondly, there is usually a difficult trade-off between 
reaching a maximum field of view and the effective pixel size 
in the resulting images. Under excellent seeing conditions, a 
critical sampling of the PSF requires pixels no bigger than 0'.'2, 
while most integral field spectrographs, including PMAS, have 
spaxels at least twice that size. As we have shown in this paper, 
image enhancement by supersampling to a finer pixel grid may 
be possible as a positive byproduct of differential atmospheric 
dispersion effects. But of course, the gain thus obtained will al- 
ways be limited by the unavoidable convolution with the true 
pixel size. 

With these provisions, a new quality level in the astrophysi- 
cal study of microlensing in QSOs can be reached. Without ma- 
jor losses in observing time, simple broad band imaging can be 
replaced by a detailed spectrophotometric investigation. In this 
paper we have presented two pieces of evidence that microlens- 
ing is important in HE 0435 — 1223: spectral differences and 
variability. Future monitoring of microlensing-induced spec- 
tral variations in this object, and in similar targets, will con- 
strain the structural properties of quasars on angular scales of 
sub-microarcseconds. Integral field spectroscopy is likely to 
become a tool of choice for this type of research. 
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